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Abstract 

Within density functional theory, we have obtained the structure of 4 He droplets doped with 
neutral calcium atoms. These results have been used, in conjunction with newly determined ab- 
initio 1 S and 1 n Ca-He pair potentials, to address the AsAp l Y\ <— 4s 2 1 So transition of the attached 
Ca atom, finding a fairly good agreement with absorption experimental data. We have studied the 
drop structure as a function of the position of the Ca atom with respect of the center of mass of 
the helium moiety. The interplay between the density oscillations arising from the helium intrinsic 
structure and the density oscillations produced by the impurity in its neighborhood plays a role in 
the determination of the equilibrium state, and hence in the solvation properties of alkaline earth 
atoms. In a case of study, the thermal motion of the impurity within the drop surface region has 
been analyzed in a semi-quantitative way. We have found that, although the atomic shift shows a 
sizeable dependence on the impurity location, the thermal effect is statistically small, contributing 
by about a 10% to the line broadening. The structure of vortices attached to the calcium atom 
has been also addressed, and its effect on the calcium absorption spectrum discussed. At variance 
with previous theoretical predictions, we conclude that spectroscopic experiments on Ca atoms 
attached to 4 He drops will be likely unable to detect the presence of quantized vortices in helium 
nanodrops. 

PACS numbers: 36.40.-c, 67.40.Yv, 33.20.Kf, 47.55.D-, 71.15.Mb 
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I. INTRODUCTION 



Optical investigations of atomic impurities in superfluid helium nanodroplets have drawn 
considerable attention in recent years.-! 2 - In particular, the shifts of the electronic transition 
lines with respect to the gas-phase transition lines (atomic shifts) are a very useful observable 
to determine the location of the foreign atom attached to the drop. Alkaline earth atoms 
appear to play a unique role in this context. While e.g., all alkali atoms reside in surface 
"dimple" states, and more attractive impurities like all noble gas atoms reside in the bulk of 
drops made of either isotope,- the absorption spectra of heavy alkaline earth atoms attached 
to 4 He drops clearly support an outside location of Ca, Sr, and Ba,-^ whereas for the lighter 
Mg atom the experimental evidence is that it resides in the bulk of the 4 He droplets.-^ 

We have recently presented Density Functional Theory (DFT) results for the structure 
and energetics of large 3 He and 4 He doped nanodroplets, showing that alkaline earth atoms 
from Mg to Ba go to the bulk of 3 He drops, whereas Ca, Sr and Ba reside in a dimple 
at the surface of 4 He drops, and Mg is in their interior." This is in agreement with the 
analysis of available experimental data, although the case of Mg has been questioned very 
recently- Moreover, according to the magnitude of the observed shifts, the dimple for alkaline 
earth atoms was thought to be more pronounced than for alkali atoms, indicating that the 
former reside deeper inside the drop than the later. This has been also confirmed by the 
calculations. In addition, the 5s5p <— 5s 2 experimental transition of Sr atoms attached 
to helium nanodroplets of either isotope has shown that strontium is solvated inside 3 He 
nanodroplets, also in agreement with the calculations.- 

Calcium atoms are barely stable on the surface of the drop, and the difference between 
the energy of the surface dimple state and that of the solvated state in the bulk of the 
drop is rather small and depends very sensitively on the Ca-He interatomic potential.— The 
aim of this work is to obtain the atomic shifts for Ca attached to large 4 He drops, and 
to compare them with the experimental data. To this end, we have improved our DFT 
approach,- treating the atomic impurity as a quantal particle instead of as an external field. 
Laser Induced Fluorescence (LIF) experiments for Ca atoms in liquid 3 He and 4 He have been 
reported 1 ^- and analyzed within a vibrating bubble model, which involves the formation of a 
bubble around the impurity, using Ca-He pair potentials based on pseudopotential SCF/CI 
calculations. 12 
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This work is organized as follows. In Sec. II we discuss the Ca-He interaction potentials 
we have used. In Sec. Ill we briefly present our density functional approach, as well as some 
illustrative results for the structure of Ca@ 4 Hejv drops. The method we have employed to 
obtain the atomic shifts is discussed in Sec. IV. In Sec. V we present the results obtained 
for calcium, discuss how thermal motion may affect the line shapes, and investigate how 
the presence of a quantized vortex line may change the Ca absorption spectrum. Finally, a 
summary is presented in Sec. VI. 



II. CALCIUM-HELIUM INTERACTION POTENTIALS 

Figure [1] shows the X X S Ca-He adiabatic potential obtained by different authors.— 
This potential determines the dimple structure described in the next Section. It can be seen 

;he others are quite similar. As in our 
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This will allow us to ascertain the 



that apart from the unpublished potential by Meyer, 
previous work,- we shall use the one obtained in Ref. 
effect of treating Ca as a quantal particle. Since the V 1 S potential seems to be fairly well 
determined, we have turned our attention to the excited adiabatic potentials. 

In a previous work 1 - the excited state potentials were calculated in a valence ab-initio 
scheme. The core electrons of calcium and helium were replaced by scalar-relativistic energy- 
consistent pseudopotentials, and the energy curves were calculated on the complete-active- 
space multiconfiguration selfconsisted field (CASSCF^^/complete-active-space multiref- 
erence second order perturbation level of theory. 

In this work we have performed fully ab-initio calculations and have only focused on 
singlet states. The calculations for excited states have been done at the CASSCF/internally 
contracted multireference configuration interaction (ICMRCI)2£>2i level of theory. In the 
calculations we have used correlation consistent polarized valence five zeta (cc-pV5Z) ba- 
sis sets. For the calcium atom we have used the (26s,18p,8d,3f,2g,lh)/[8s,7p,5d,3f,2g,lh] 
basis set developed by Koput and Peterson,— and for the helium atom we have used the 
(8s,4p,3d,2f,lg)/[5s,4p,3d,2f,lg] basis set developed by Woon and Dunning.— 

The calculations were performed within the MOLPRO suite of ab-initio programs. 24 
The molecular orbitals used for the excited states calculations were optimized in the state 
averaged CASSCF method for all singlet states correlating to (4s 2 ) 1 S, (4s3d) 1 D, and (4s4p) 1 P 
atomic asymptotes. The active space was formed by distributing the two valence electrons 
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of the Ca atom into As3d4p valence orbitals. The Is2s2p3s3p orbitals of calcium and Is 
orbital of helium were kept doubly occupied in all configuration state functions, but they 
were optimized in the CASSCF calculations. The resulting wave functions were used as 
references in the following ICMRCI calculations. At the ICMRCI level, the Is2s2p3s3p 
orbitals of calcium and Is orbital of helium were kept as doubly occupied in all reference 
configuration state functions, but these orbitals were correlated through single and double 
excitations. The 4s, 4p, and 3d calcium orbitals had not restricted occupation patterns. We 
show in Fig. [2] the excited adiabatic potentials; to have better insight into the potential 
minima of the X S and 1 I1 potentials, they have been plotted correlating to the (4p) x P Ca 
term. 



III. DFT DESCRIPTION OF HELIUM NANODROPLETS 

In recent years, static and time-dependent density functional methods 2 ^^ have become 
increasingly popular to study inhomogeneous liquid helium systems because they provide 
an excellent compromise between accuracy and computational effort, allowing to address 
problems inaccessible to more fundamental approaches, see e.g. Ref. y for a recent review. 
Obviously, DFT cannot take into account the atomic, discrete nature of these systems, but 
can address inhomogeneous helium systems at the nanoscale^S and take into account the 
anisotropic deformations induced by some dopants in helium drops. Both properties are 
essential to properly describe these systems. 

Our starting point is the Orsay-Trento density functional,— together with the Ca-He 



adiabatic potential X X S of Ref. [15|, here denoted as Vca-He- This allows us to write the 
energy of the Ca-drop system as a functional of the Ca wave function $(r) and the 4 He 
"order parameter" \l/(r): 



= — — / d 3 r|W(r)| 2 + / d 3 r£(p) 
2 rriHe J J 

+ ^— [ d 3 r|V$(r)| 2 + / / d 3 rdV|$(r)| 2 ^ Me (|r-r'|)p(r') . (1) 
2mca J J J 

The order parameter is defined as \l/(r) = J p(r) exp[«S'(r)], where p(r) is the particle density 

and v(r) = ^VS , (r)/m 4 is the velocity field of the superfluid. In Eq. ([1]), £(p) is the 4 He 

"potential energy density" .— In the absence of vortex lines, we set 5* to zero and E becomes 
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a functional of p and $. Otherwise, we have used the complex order parameter ^(r) to 
describe the superffuid. 

We have solved the Euler-Lagrange equations which result from the variations with re- 
spect to ty* and $* of the energy E[$f, <&] under the constrain of a given number of helium 
atoms in the drop, and a normalized Ca wave function, namely: 

'"' A^ + \ 6 -f + U He )y = ^ (2) 



2m He [5p 

h 2 



A<5> + U Ca $> = e<S> , (3) 
2m Ca 

where p is the helium chemical potential and e is the lowest eigenvalue of the Schrodinger 
equation obeyed by the Ca atom. The effective potentials Uue and Uc a are defined as 



U He (r) = |dV|$(r')| 2 Vc a „ He (|r-r'|) 

U Ca (v) = J dVp(r') V Ca . He {\v - r'|) . (4) 

The coupled Eqs. (J2H3J) have to be solved selfconsistently, starting from an arbitrary but 
reasonable choice of the unknown functions ^ and $. In spite of the axial symmetry of the 
problem, we have solved them in three-dimensional (3D) cartesian coordinates. The main 
reason is that these coordinates allow us to use fast Fourier transformation techniques 29 to 
efficiently compute the convolution integrals entering the definition of S(p), i.e. the mean 
field helium potential and the coarse-grained density needed to compute the correlation term 
in the He density functional,— as well as the fields defined in Eqs. (jl]). 

The differential operators in Eqs. (JT}{3j) have been discretized using 13-point formulas for 
the derivatives, and Eqs. (GHH]) have been solved employing an imaginary time method; 30 
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Typical calculations have been 



some technical details of our procedure are given in Ref. 
performed using a spatial mesh step of 0.5 A. We have checked the stability of the solutions 
against reasonable changes in the step. 

Equations ([2113]) have been solved for several N values from 100 to 2000. They will allow 
us to study the atomic shift as a function of the cluster size. The equilibrium configurations 
of Ca@ 4 Heiooo and Ca@ 4 He2ooo will be shown later on. 

Figure [3] shows the energy of a Ca atom attached to a drop, defined as the energy difference 
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S N (Ca) = E(Ca@ 4 He N ) - E( 4 He N ) . (5) 

On the figure are shown also the results obtained treating calcium as an external field. 8 It 
can be seen that for large drops the energy of the calcium atom is about 10 K less negative 
due to its zero point motion. We want to stress again how barely stable is the calcium atom 
on the surface of 4 Hejv drops. For instance, we have found that the total energy of the 
equilibrium -dimple- configuration of Ca@ 4 Hei oo is ~ —5467.4 K, whereas it is ~ —5455.0 
K when Ca is forced to be at the center of the drop. For Ca@ 4 He5oo, the corresponding 
values are ~ —2525.2 K and ~ —2511.1 K, respectively. 

The dimple depth £, defined as the difference between the position of the dividing surface 
at p = po/2, where po = 0.0218 A -3 is the bulk liquid density, with and without impurity, 
respectively, is shown in Fig. HI Due to the zero point motion that pushes the impurity 
towards lower helium densities, for large drops the dimple depth is about 0.8 A smaller when 
the zero point motion is included than when it is not. 8 This change in the depth is large 
enough to produce observable effects in the calculated absorption spectrum, as discussed 
below. 

It can be seen that the dimple depth curve £(iV) has some structure. This is not a 
numerical artifact, but a genuine effect due to the interplay between the Ca atom and the 
drop, whose density, even for pure drops, shows conspicuous oscillations all over the drop 
volume, extending up to the surface region irrespective of whether the drop is described 
within DFT or Diffusion Monte Carlo met hods .-^'22' The interplay of these oscillations with 
those arising from the presence of the impurity little affects the total energy of the system 
and hence, the Ca energy, but yields some visible structure in the density distributions that 
shows up in related quantities, like the dimple depth. To illustrate it, we show in Fig. [5] the 
density of the helium moiety of Ca@ 4 He2oocb where the interference pattern can be clearly 
seen. 

Further insight can be gained studying, for a given drop, the energy of the Ca@ 4 He7v 
complex as a function of the distance between the centers of mass of the impurity and of 
the helium moiety. This can be done adding an appropriate constraint to the total energy 
in Eq. ([1]), and solving the corresponding Euler-Lagrange equations. Specifically, we have 
minimized the expression 
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E + f [Z-Zof 



(6) 



where Z is the average distance in the z direction between the impurity and the geometrical 
center of the helium moiety 



and Xc is an arbitrary constant, large enough to guarantee that upon minimization, Z equals 
the desired Z value. We have also applied this method to Mg doped helium drops, and will 
present the details of the calculation elsewhere.— 

We show in the bottom panel of Fig. [6] the total energy of the Ca@ 4 Hesoo and Ca@ 4 Heiooo 
systems as a function of Zq. For the sake of comparison, the energies and Zq distances are 
referred to their equilibrium values. The vertical lines roughly delimit the drop surface 
regions, conventionally defined as the radial distance between the points where the density 
equals O.lpo and 0.9po (see Fig. [TUI) . The horizontal line has been drawn 0.4 K above the 
equilibrium energy, representing the accessible energy range due to the temperature of the 
helium drops. 34 Its intersection with the energy curves yields a qualitative measure of the 
dispersion of the impurity location due to thermal motion. The energy curve of Ca@ 4 Heiooo 
displays some structure to the left of the minimum due to the mentioned interference pattern. 
This behavior has not been disclosed before, and appears in the course of "pushing" the 
impurity inside the droplet from its equilibrium position (a similar structure shows up for 
Ca@ 4 He5oo, but at more negative A^o values). While it affects rather little the equilibrium 
location of the Ca atom because of its clear surface location, and hence the atomic shift -see 
however the inset in Fig. [8j it plays a substantial role in the solvation of magnesium atoms 
in small helium drops. 33 We recall that this problem has recently drawn the attention of 
experimentalists and theoreticians as well.— iZ&MiL2& How the position of the Ca atom affects 
the absorption spectrum throughout the change in the dimple structure will be discussed in 
Sec. V. 

IV. EXCITATION SPECTRUM OF AN ATOMIC IMPURITY IN A 4 HE DROP 

Lax method 37 offers a realistic way to study the absorption spectrum of a foreign 
atom embedded in helium drops. It makes use of the Franck-Condon principle within a semi- 

7 
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classical approach, and it has been employed to study the absorption spectrum of several 
atomic dopants attached to fairly small 4 He drops.— The case of alkali atoms attached 
to large drops described within DFT has been also considered, see Refs. I4l|j42l and refer- 
ences therein. Lax method is usually applied in conjunction with the diatomics- in- molecules 
theory— which means that the atom-drop complex is treated as a diatomic molecule, where 
the helium moiety plays a role of the other atom. 

In the original formalism, to obtain one has to carry out an average on the possible 
initial states of the system that may be thermally populated. Usually, this average is not 
needed for helium drops, as their temperature, about 0.4 K,— is much smaller than the 
vibrational excitation energies of the Ca atom in the mean field represented by the second 
of Eqs. (jlj) . 44 However, thermal broadering due to the "wandering" of the dopant must be 
analyzed separately if it is in a dimple state, as this may have some influence on the line 
shape. Its effect on the absorption spectrum will be exemplified later on for the N = 500 
case. Contrarily, when the impurity is very attractive and resides in the bulk of the drop, 
thermal motion plays no role, as the dopant hardly gets close enough to the drop surface 
to have some effect on the line shape. In this case, dynamical deformations of the cavity 
around the impurity may be relevant (Jahn- Teller effect). It cannot be discarded that, if 
some of these very attractive impurities have an angular momentum large enough,- 4 ^ they 
may get close to the drop surface, in which case thermal motion might have some influence 
on the absorption spectrum. 

We review here the essentials of the method and the way we have implemented it. In 
particular, we present some of the expressions in cartesian coordinates, better adapted to 
our approach. They are of course equivalent to the expressions in spherical coordinates that 



can be found in the literature -see e.g. Ref. 



39 and references therein. 



A. Line shapes 

The line shapes for electronic transitions from the ground state (gs) to the excited state 
(ex) in a condensed phase system can be written as 

I(u) oc j dt e- <wt (* fl "|Dj e et«- D ge e -^'|* flS ) , (8) 
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where D ge is the matrix element of the electric dipole operator, H. gs and 7i ex are the Hamil- 
tonians which describe the ground and excited states of the system respectively, and \ty gs ) 
represents the ground state. can be evaluated using the Born-Oppenheimer approxima- 
tion, which makes a separation of the electronic and nuclear wave functions = \e l )\ip l ), 
and the Franck-Condon principle, whereby the heavy nuclei do not change their positions 
or momenta during the electronic transition. If the excited electron belongs to the impurity, 
the helium cluster remains frozen, so that the relevant coordinate is the relative position 
r between the cluster and the impurity. That principle amounts to assuming that D ge is 
independent of the nuclear coordinates. Taking into account that e~^ H9a \^ 9S ) = e~ ltul3a \^ 9S ) 
and projecting on eigenstates of the orbital angular momentum of the excited electron |m), 



I(u) oc E / dt e- l{uj+LVgs)t (^ es \Dl\m) e~% n ™ (m\D gc \q ss ) 

= E I dt e-^ +£jgs )' f d 3 r f d 3 r' {^ s \r){r\D\ c \m)e^ {m\D gc \v'){r'\^ s ) 
m J J J 

= \D gc \ 2 E / dt e- i(tJ+ ^ S) * / d 3 r (r)* e* *™ (r) , (9) 

where hu>x and ipx ( r ) are the energy and the wave function of the ro- vibrational ground state 
of the frozen helium-impurity system, and H^(r) is the ro-vibrational excited Hamiltonian 
with potential energy V^ x (r) determined by the electronic energy eigenvalue, as obtained 
in the next subsection for a p <— s transition. Eq. (Q will be referred to as the Fourier 
Formula, and it is the Fourier transform of the time-correlation function. It is nothing but 
a sum at the resonant energies weighted with the well-known Franck-Condon factors: 



I(u) oc E l dt e-^-^f I |Vf ) 

m 

= EE / dt ^^(^mm mmi$) 

m v u' 

= EE/ d ^^ f ^KCI^f)l 2 

m v 

= EE^ + ^f-OKCI^f)! 2 , 

m v 

where cu™ and are the ro-vibrational eigenvalues and eigenstates of the Hamiltonian 

rrcx 
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If the relevant excited states for the transition have large quantum numbers, they can be 
treated as approximately classical— using the averaged energy hu™ ~ V^(r) which is 
independent of v. In this case we obtain the expression 

cx W d 3 r |Vf (r)| 2 % + u,f -V™{r)/K) 

m 

= K f cfr '^' r " 2 (10) 

where VL m {uj) is the surface defined by the equation to + Ux — Krf( r )/^ = 0- We will refer 
to this equation as the Semi- Classical Formula. 

If the atom is in bulk liquid helium, or at the center of the drop, the problem has spherical 
symmetry and the above equation reduces to 

oc 4tt]T I dr \r ^(r)\ 2 5(uj + uj g x - V^(r)/h) 
[r (r)] 2 



dV m ex (r)/dr 



(11) 

r=r m (a;) 



where r m {uj) is the root of the equation uj + Ux — Kn X ( r )/^ = 0- 

In the non-spherical case, we have evaluated from the first expression in Eq. (TTU|) 
using the discretization 



J(nAw + w ) oc YlY^\^x (A*ijk + *o] 



2 AxAyAz 



x {6 [{n - 1/2) Alu + uj + w f - V; e t x (Ar ijfe + r )/ft] 

[(n + 1/2) Acu + cu + - ^(Ar^fc + r )/ft]} , (12) 

where G is the step function and Auj is a frequency step small enough so that the above 
discretization represents the delta function. We also take advantage that only points near 
the impurity contribute to the integral, by writing Ar^+ro = (iAx+x , jAy+y , kAz+z ), 
with Ax, Ay, Az being the spatial mesh steps used in the discretization, and being r an 
arbitrary point in the neighborhood of the impurity. Finally, we recall that I(u>) needs to be 
evaluated only in a narrow frequency range starting from an arbitrary Uq which can be, e.g., 
the free atom frequency. This range defines the maximum n value in the above equation. 
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B. Excited ro-vibrational potential for a p <— s transition 

We next determine the potential energy surfaces V^(r) needed to carry out the calcula- 
tion of the atomic shifts. 

1. Pairwise sum aproximation 

The pair-interaction between an atom in a s-state and an atom in a p-state can be 
expressed, in the cartesian eigenbasis (\x), \y), \z)) as 



( V u (r) N 



U{r) 



V u (r) 
V Vs(r) J 
= V n {r){\x){x\ + \y){y\} + Vv(r)\z){z\ 

= Vn(r)I + {Vx(r)-V n (r)}\z){z\ , (13) 

where Vn{r) and Vs(r) are the adiabatic potentials neglecting the spin-orbit interaction and 
r is the distance between atoms. For a system of N helium atoms and an excited impurity 
in a p-state, the total potential is approximated by the pairwise sum 



^=E{^n(r„)I+[\/ s (r„)-Vri(r ri )] J R„|^)(^| J R; 1 } , (14) 

71=1 

where r„ is the distance between the n th helium atom and the impurity, and R n is the 
rotation matrix that transform the unity vector z into the f n vector. It can be shown that, 
in cartesian coordinates, 

(*VU*><*KV> = ^, (is) 

where x 1 = x, x 2 = y, x 3 = z, and r 2 n = x 2 n + y 2 + z 2 . Thus, the matrix elements of the total 
potential are 

N ( i j ) 

( x *\U\xi) = U tj = E VnirJSij + [V^r n ) - V u (r n )}^^ . (16) 

n=l I r n ) 

Using the continuous density approach inherent to DFT [J2 n ~^ I d 3 r'p(r')], this expression 
can be written as 
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UM = J dVp(r' + r) |Vn(r')% + [Vs(r') - ^n(r')]^| • (17) 

The eigenvalues of this symmetric matrix are the sought-after V^ x (r) which define the po- 
tential energy surfaces (PES) as a function of the distance between the centers of mass of 
the droplet and of the impurity, and are given by the three real roots Aj(r) of the equation 



A 3 + C\ 2 + BX + A = (18) 



with 



c = -J2u u 



2 

2 z < \ " I J ~ ij 

3/1 ] •) 



A = E {^UiiUf k --U ii U jj U kk --U ij U jk U ki ) . (19) 

It can be shown that for spherical geometry, Eq. (TP7l) is diagonal with matrix elements 
(in spherical coordinates) 



Xi(r) = Uii(r) = 2tt J J r' 2 sin 6'd6'dr' p(\r' + r + 2rVcos0'|) 

"I 



x \Vn(r') + lVx(r')-V u (r')} 



-{Sn + 5 i2 )sin 2 e' + S i3 cos 2 & 



(20) 



2. Spin- Orbit coupling 



For atomic impurities in which the spin-orbit (SO) interaction is prominent and compa- 
rable to the splitting of the P-states due to the interaction with the droplet, it has to be 
taken into account in the calculation of the PESs. This is usually done considering that 
the SO splitting of the dopant is that of the isolated atom irrespective of the impurity-drop 
distance.— Given the atomic structure of alkaline earth atoms, the SO interaction can 
be safely neglected in the PES calculation. However, we discuss it here for the sake of 
completeness and future reference. 

When the spin-orbit interaction is taken into account, the total potential can be 
written as = U + Vso, where V$o has the form, in the spin-cartesian orbit basis 
(\x, 1/2), \x, -1/2), \y, 1/2), \y, -1/2), \z, 1/2), \z, -1/2)): 
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Vf 



A, 



so 



-% 1 

i -1 

% -i 

-i -i 

-1 i 

10 i 



where A^ s is 2/3 of the experimental SO splitting of an isolated atom. Kramers' theorem 
states that there is a two-fold degenerate manyfold of systems with a total half-integer spin 
value that cannot be broken by electrostatic interactions,— so that the two-fold degenerate 
eigenvalues that define the PES's are the roots of the equation 



A 3 + C\ 2 +[B- ltf u ]\ + {A + \[Al - AlC}} = (21) 

with A, B and C defined in Eq. (fl9l) . In the case of spherical geometry, the eigenvalues 
adopt a simple expression: 



Ai(r) 

A 2 (r) 
As(r) 



\(U 11 + U 33 ) + ± 



-A ls + J9AI - 4A £s (Un - U 33 ) + 4{U U 



u 33 y 



u 



ii 



2 



2(^11 + ^33) 



A es - J9AI - 4A es (U u - U 33 ) + 4(U U - U 33 y 



(22) 



which reduces to Eq. f l20l if A ia = 0. 

It is important to notice that for spherical geometries (spherically symmetric impurities 
in liquid helium or at the center of a drop), when the SO interaction is negligible two of the 
PES are degenerate, as it can be seen from Eq. ( 1201) . that yields Ai(r) = A2(r) 7^ As(r).— 
Thus, the existence of the SO interaction not only is the reason of the appearance of the D\ 
and D 2 lines in the case, e.g., of alkali atoms in bulk liquid helium, but also the reason of 
either the broadening or the splitting of the D 2 line.— For this particular geometry, another 
contribution to the splitting of the D 2 line is the Jahn- Teller effect caused by dynamical 
quadrupole deformations of the cavity surrounding the impurity, which develop irrespective 
of whether the spin-orbit energy is relevant or not.- 1 ^ When the impurity resides in a 
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deformed environment like a dimple, the three PES are non degenerate, and may cause 
the appearance of three distinct peaks in the absorption spectrum, or of just one single 
broad peak, as it happens in the case of Ca, Sr and Ba atoms attached to 4 He drops.- 1 - The 



liquid 4 He results for alkaline earth atoms are reported in Refs. 



52 



53 and references therein. 



V. RESULTS FOR THE ABSORPTION SPECTRUM OF CALCIUM ATOMS 



A. Line shifts 



The problem of obtaining the line shifts has been thus reduced to that of the dopant in 
the 3D trapping potentials corresponding to the ground state, Uc a ( r ), an d P excited states, 
Aj(r). Since we have neglected the fluctuations of the dimple -shape fluctuations^ 4 - and their 
coupling to the dopant dipole oscillations, as well as inhomogeneous broadening resulting 
from droplet size distributions, laser line width and similar effects, the model is not expected 
to yield the line shapes, but only the energies of the atomic transitions. These limitations 
are often overcome by introducing line shape functions or convoluting the calculated lines 
with some effective line profiles.— 1 ^ We discuss now some illustrative examples without 
considering these justified but somewhat uncontrolled convolutions. 

Figure [7] shows the absorption spectrum of Ca@ 4 He2ooo calculated with the the Semi- 
Classical Formula. The much involved Fourier Formula calculation is unnecessary in the 
Ca@ 4 Hejv case. The reason is the absence of bound-bound transitions from the ground 
state PES to the II or E ones because their wells are spatially well apart. The starred 
vertical line represents the gas-phase transition. The three components of the absorption 
line, each arising from a different excited PES, are also shown. We have normalized to one 
the integral of each component. This choice comes out naturally from the normalization of 
the wave function of the impurity; obviously, the relative intensity of the three components 
is not arbitrary. In Fig. [7J the n PESs contribute to build up the maximum of the line, 
whereas the the long blueshift tail arises from the S PES. 

No appreciable differences appear between N = 2000 and 2500, the largest drop we have 
calculated. This saturation has been also observed in the experiment,-^ although for larger 
mean cluster sizes, about iV ~ 3000. The peak energy is 79 cm -1 , a 10% larger than the 
experimental saturation value of 72 cm" 1 ,- which indicates a fairly good agreement between 
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theory and experiment. The calculated width (FWHM) is ~ 55 cm -1 , whereas the measured 
width for drops in the N = 1500 — 2000 range is ~ 140 — 155 cm" 1 ,- i.e., about three times 
larger. We will show later on how thermal motion affects the theoretical result. 

Figure [S] shows the total absorption spectrum of Ca attached to 4 Hejv droplets for several 
N values. It is interesting to notice the evolution of the absorption line as the number of 
atoms increases. As a general rule, the smaller the drop, the smaller the splitting of the II 
components. Indeed, they would be degenerate if N = 1, as Eq. ([TBI shows. This explains 
why the main peak for N = 100 is the narrower one, and is the reason why the main peak 
is fairly apart from the £ shoulder. As N increases, so it does the splitting, while the three 
components of the peak become broader. Eventually, if N is large enough, it is not possible 
to distinguish the components of the absorption line. It is also obvious that line broadening 
due to effects not considered here may wash out the blueshifted shoulder found for small N 
values. 

The inset in Fig. [8] shows the calculated shifts relative to the gas-phase transition, 
compared with the experimental values.- One can appreciate a small oscillation for the 
largest drops; it is a genuine effect produced by the dimple structure. Further insight can be 
gained from the study of the absorption spectrum as a function of Z . To this end, we display 
in the top panel of Fig. [9] the absorption spectra for the N = 500 droplet corresponding to 
Z values from 16 to 19 A in 0.5 A steps, all of them within the drop surface region (the 
equilibrium value is Zq= 17.45 A). The inset in Fig. [9] shows the dependence of the relative 
shift on the location of the Ca atom. One can see how a 3 A dispersion in the position of 
the impurity generates a ~ 45 cm -1 change in the shift, showing in a quantitative way the 
well known sensitivity of this quantity to the structure and depth of the dimple. It is worth 
seeing how the shift decreases as the distance between the centers of mass increases, and at 
the same time the absorption peak becomes more asymmetric as the Ca environment does 
(see also Subsection C). 

B. Thermal broadening 

As we have indicated, for very attractive impurities, the foreign atom is fully solvated, 
and its derealization in the bulk of the drop due to thermal motion hardly introduces a 
significant change in the line shape. When the dopant is in a dimple state, it has to be 
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checked whether thermal motion may have observable effects on the absorption spectrum, 
as Fig. [6] seems to indicate for calcium. 

To ascertain this effect, we have carried out a thermal average of the spectrum using an 
approximate expression for the probability density. Referring the energy Ej, of a given Z Q . 
configuration to the equilibrium value, AEi = Ei — E gs , and neglecting the kinetic energy 
of the impurity and the displaced fluid, we write the probability density for the position Z^ i 
of the Ca atom as 

^2 e -AEi/k B T 

Wl = E 3 -Z2.e-^./fcBT AZo . ' ( 23 ) 
where ks is the Boltzmann constant, T = 0.4 K, and the sum -actually integral- runs on 
the selected Z . configurations. 5 - The Zq. factor takes care of the relative volume available 
to each configuration. With this definition the probability of finding the Ca atom between 
Zo i and Zq { + AZo t is WiAZ^. We show in the top panel of Fig. [6] the probability densities 
Wi corresponding to the configurations displayed in the bottom panel. We see that there is 
a non-negligible probability of finding the Ca atom in a broad region of the drop surface, 
and consequently we have addressed, of study, the statistical properties of the 

Ca@ 4 He 50 o system at this temperature. 

We have found that the mean position, calculated as (Z ) = I±iZ 0i WiAZ Oi , and the stan- 
dard deviation, calculated as u(Z ) = yj (Zq) — (Z ) 2 are 17.38 A and 0.76 A, respectively. 
This dispersion in the position generates a dispersion in the value of the shift. To quantify 
this effect we have evaluated the mean value of the shift, calculated as (Au) = Y^AuJiWiAZQ., 
and its standard deviation, calculated as ct(Alu) = \J (Alu 2 ) — (Alu) 2 , using the values shown 
in the inset of the top panel of Fig. [9j which correspond to the seven lower energy configu- 
rations in Fig. [6l We have obtained (Au) ± a(Au) = 63.8 ± 11.5 cm"" 1 . 

The thermally averaged Ca absorption spectrum is shown in the bottom panel of Fig. [9] 
(solid line), as well as that corresponding to the equilibrium configuration (dashed line). To 
carry out the average, we have used the in the top panel of Fig. [91 and averaged them 
as I(oS) = T,iIi(uj)wiAZ 0i . This procedure, consistent with the Franck-Condon principle, as- 
sumes that absorption proceeds instantaneously on any of the frozen drop-Ca configurations 
characterized by a Z^ i value. 

We are led to conclude that the thermal motion effect is rather small. It increases the 
FWHM by about 10%, from ~ 49.5 cm" 1 to ~ 55.0 cm -1 , still a factor of three smaller than 
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the experimental value. From Fig. [6j we expect a similar effect for the N = 1000 drop, and 
likely for larger drops. 



C. Calcium atoms attached to vortex lines in 4 He drops 

Since 4 He is superfluid, it is quite natural to wonder about the appearance and detection 

571 and references therein. Adapting an 



of quantized vortices in droplets, see e.g. Refs. 
idea originally put forward by Close et al.,— it has been proposed^ that Ca atoms should 
be the dopant of choice to detect vortices by means of microwave spectroscopy experiments. 
The rationale of this proposal is that Ca atoms are barely stable on the drop surface and 
become solvated in its interior in the presence of a vortex line.— These conclusions were 
drawn from DFT calculations using Meyer's Ca-He potential which, as shown in Fig. [TJ is 
slightly more attractive than recent potentials. If this scenario were plausible, one would 
not need the microwave spectroscopy experiments suggested in Ref. to detect a vortex 
state in a Ca@ 4 Hejv drop: LIF spectroscopy could do the job, given the sizeable difference 
between the blueshifts of the absorption lines when Ca has been drawn inside the drop by 
the vortex (similar in value to the liquid helium blueshift), and when it resides in a dimple 
state in vortex-free drops.- 

This has prompted us to re-analyze the structure of a large N = 1000 drop hosting a 
calcium atom attached to a vortex line along the symmetry axis. Within DFT, a robust 



method to generate vortex configurations in liquid helium is described in Ref. [59|. We adapt 
it here to the case of helium drops. For a n — 1 quantum circulation vortex line, we start 
the imaginary time evolution to solve Eq. @ from the initial state 

P 1/2 (r) 

^(r) = ^=LL( X + l y) (24) 

if x and y are nonvanishing, and zero otherwise, where p(r) is the vortex-free Ca@ 4 Heiooo 
helium density. After the minimization procedure is converged, we have checked that the 
obtained final configuration is indeed a n = 1 vortex state. 

Figure [10] shows equidensity lines for the equilibrium configuration of Ca@ 4 Hei oo with 
and without a vortex line along its symmetry axis. It can be seen that the vortex line draws 
the impurity towards the bulk of the droplet, but it still resides in a deeper surface dimple. 

Figure [H] shows the absorption spectrum for the two configurations displayed in Fig. [TUl 
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The effect of the presence of the vortex on the absorption spectrum of calcium is twofold. On 
the one hand, the maximum of the absorption peak is shifted towards the bulk value because 
of the deeper dimple. On the other hand, the FWHM increases by about a factor of two. 
The reason is the spreading of the Ca wave function within the stretched Uc a { r ) well, that 
allows the atom to "probe" a wider region in the excited PES, thus increasing the width. 
Notice also the larger splitting of the peaks that form the line due to the more anisotropic 
helium environment. This also contributes to increasing the width of the absorption peak. 
Unfortunately, the experimental absorption line is so broad and asymmetric that the extra 
shift caused by the vortex is not enough to displace the line to a region where it could be 
distinguishable on top of the vortex-free absorption line. 

VI. SUMMARY 

Within density functional theory, we have carried out a detailed study of the absorption 
spectrum of calcium atoms attached to 4 HeAr drops in the vicinity of the AsAp 1 Pi <— 4s 2 
1 So transition, finding a semi-quantitative agreement with experiment. To this end, we have 
improved our previous implementation of the DF method by incorporating the zero point 
motion of the impurity, and have carried out ab-initio calculations to obtain the excited X S 
and 1 I1 Ca-He potentials needed to obtain the potential energy surfaces. 

We have studied the drop structure, finding that the "interference" between the density 
oscillations of the helium moiety arising from its intrinsic structure and those arising from 
the presence of the impurity plays a role in the determination of the position of the impurity. 
This may be relevant for the solvation of alkaline earth atoms, especially for magnesium.— 

In a case of study, we have systematically addressed the dependence of the relative shift 
on the position of the impurity, quantitatively assessing the relevance of a proper description 
of the dimple to reproduce the experimental results. We have statistically taken into account 
the influence of the thermal motion of the impurity on the absorption line, concluding that 
it only increases the line width by a modest amount. 

Finally, we have addressed the Ca absorption spectrum when the helium drop hosts a 
vortex line, and conclude that absorption spectroscopy experiments on these drops would be 
likely unable to ascertain the presence of vortical states. In spite of this, the study of vortex 
lines pinned by calcium atoms in superfluid helium drops is interesting by itself, especially 
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the evaluation of the atomic shift caused by the presence of a vortex line. 
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FIG. 1: (Color online) Some recent X X S Ca-He pair potentials: (1) Ref. [l3|; (2) Ref. 
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FIG. 2: (Color online) Excited Ca-He pair potentials used in this work 
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FIG. 3: (Color online) Energy (K) of a calcium atom as a function of the number of 4 He atoms in the 
droplet (dots). Results obtained treating calcium as an external field are also shown (triangles). 8 
The lines have been drawn to guide the eye. 
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FIG. 4: (Color online) Calcium dimples (A) as a function of the number of 4 He atoms in the 
droplet (dots). Results obtained treating calcium as an external field^ are also shown (triangles). 
The lines have been drawn to guide the eye. 
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FIG. 5: (Color online) Helium density of the Ca@ 4 He2ooo drop on the x = plane. The probability 
density of the calcium atom is also shown, rescaled multiplying it by a pq factor for the sake of 
clarity. 
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FIG. 6: (Color online) Bottom panel: total energy (K) of Ca@ 4 Hesoo (circles) and Ca@ 4 Heiooo 
(squares) as a function of Z$ (A) , both referred to their equilibrium values. The vertical lines delimit 
the drop surface regions, and the horizontal line has been drawn 0.4 K above the equilibrium 
energy. Top panel: probability densities for the configurations displayed in the bottom panel. 
Circles correspond to Ca@ 4 Hesoo, and squares to Ca@ 4 Heiooo- in both panels, dotted lines refer 
to the N = 500 drop, and dashed lines to the N = 1000 drop; the lines (cubic splines) have been 
drawn to guide the eye. 
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FIG. 7: (Color online) Calcium absorption spectrum for N = 2000 in the vicinity of the 4s4p 
1 Pi <— 4s 2 1 So transition. The starred vertical line represents the gas-phase transition, and the 
dotted vertical line the experimental value extracted from Fig. 3 of Ref. [4j. The three components 
of the absorption line, each arising from a different excited PES, are also shown. 
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FIG. 8: (Color online) Total absorption spectrum of Ca attached to 4 Hejv droplets in the vicinity of 
the 4s4p 1 Pi <— 4s 2 1 So transition, for the indicated N values. The starred vertical line represents 
the gas-phase transition. The inset shows the calculated shifts relative to the gas-phase transition 
(dots), defined as the energy of the maximum of the absorption line minus the energy of the gas- 
phase transition. The experimental values are also shown (triangles). The arrow indicates the 
asymptotic experimental value.— 
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FIG. 9: (Color online) Bottom panel: Thermally averaged absorption spectrum of Ca in a N = 
500 drop (solid line); the dashed line represents the absorption spectrum corresponding to the 
equilibrium configuration. Top panel: Ii(to) spectra used to carry out the average. The Zo i values 
go from 16 to 19 A in 0.5 A steps. The inset shows the associated atomic shifts. 
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FIG. 10: (Color online) Equidensity lines showing the equilibrium configuration of a Ca atom on 
a 4 Heiooo droplet with a vortex line along its symmetry axis (top panel), and without it (bottom 
panel). The lines correspond to densities 0.9po to O.lpo m O.lpo steps; these lines span the drop 
surface region as well as the vortex core. The equidensity lines of the Ca probability density are 
similarly plotted starting from its maximum value. 
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FIG. 11: (Color online) Ca absorption spectrum for the N = 1000 droplet of Fig. [10] with (dashed 
lines) and without a vortex line along its symmetry axis (solid lines). The absorption line has 
been decomposed into its three components. The starred vertical line represents the gas-phase 
transition, and the dotted vertical line represents the experimental value for bulk liquid 4 He?ii 
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